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Abstract 



Oh! 

^ • We combine the constraint suggested by the recent BNL E821 measurement of the 

anomalous magnetic moment of the muon on the parameter space of the constrained MSSM 
^ I (CMSSM) with those provided previously by LEP, the measured rate of 6 57 decay and 

the cosmological relic density Q^h'^. Our treatment of f^^^^ includes carefully the direct- 
channel Higgs poles in annihilation of pairs of neutralinos x a complete analysis of x — ^ 
coannihilation. We find excellent consistency between all the constraints for tan /5 > 10 
and /i > 0, for restricted ranges of the CMSSM parameters mo and mi/2. All the preferred 
CMSSM parameter space is within reach of the LHC, but may not be accessible to the Teva- 
tron collider, or to a first-generation e~^e~ linear collider with centre-of-mass energy below 
1.2 TeV. 
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The recent BNL E821 measurement [0] of the anomalous magnetic moment of the muon, 
a,, = {g^ — 2)/2, may indeed be a harbinger of new physics [0] beyond the Standard Modeh 



6a 



a 



exp 



a 



SM 



(43 ± 16) X 10^^*^. The largest error in the Standard Model prediction 
is that due to the hadronic contributions, principally vacuum polarization diagrams, with 
the most important uncertainty being that in the low-energy region around the p peak. 
The value of these hadronic contributions ^ used in the E821 paper IQ does not include 
the latest data from Novosibirsk 0], Beijing and CLEO 0, but these are unlikely to 
change the overall picture: we recall that the hadronic error < 7 x 10^^° is much smaller than 
the apparent discrepancy and the experimental error. Advocates of new physics beyond the 
Standard Model may therefore be encouraged. However, we recall that the Z ^ bb branching 
ratio was once thought to show a bigger discrepancy with the Standard Model, and we also 
caution that the 2.6 cr significance of the muon anomaly is formally less than the preliminary 
2.9 cr significance of the LEP Higgs 'signal' [^. 

A priori, the BNL measurement favours new physics at the TeV scale, and we consider 
the best motivated candidate to be supersymmetry. Even before the hierarchy motivation for 
supersymmetry emerged, the potential interest of was mentioned, and a pilot calculation 
performed Soon after the realization that supersymmetry could alleviate the hierarchy 
problem, the first 'modern' calculations of supersymmetric contributions to were pub- 
lished |TD|, [IT], [T2|, |T^. These were followed by more complete calculations |T^, |T3|, [T^, 
including the mixing expected for neutralinos, charginos and smuons. In particular, it was 
noted in |T7] that some contributions are enhanced at large tan/3. The supersymmetric 
calculations we use in this paper are taken from - for other recent calculations, see [[T9| , 



and we include the leading two-loop electroweak correction factor pOf . For some time, it has 
been emphasized [|l7i [2T|] that the BNL experiment would be sensitive to a large range of 
the parameter space of the constrained minimal supersymmetric extension of the Standard 
Model (CMSSM) with universal soft superymmetry-breaking parameters at the input GUT 
scale, determining in particular the sign of the Higgs mixing parameter fi |1T7|, |2l||. Combining 
these calculations with the BNL measurement, /z > is favoured, along with values of tan/5 
that are not very small. 

The constraints from the E821 experiment are particularly interesting when combined 
with the information from LEP P], the measured value of the b ^ s'y decay rate and 
restrictions on cold dark matter imposed by astrophysics and cosmology, assuming that 



the lightest supersymmetric particle (LSP) is the lightest neutralino x ll^l) ^ind that R 
parity is conserved. Several combinations of these other constraints have been made by 
us 0, ^, H, 13 0, g and others [| 



31], before the advent of the E821 result. 
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We draw particular attention to a recent combined analysis p8[ of these constraints 



at large tan/3 > 20, which benefited from recently available b ^ s^y calculations and 
made new calculations at large tan/5 of the the relic density Q^h"^. We found |^ two 
important effects on the calculation of Q^h'^, due to improvements of previous calculations of 
coannihilations and direct-channel XX annihilations through the heavier neutral MSSM 
Higgs bosons H and A. Both of these effects extended the region of CMSSM parameter space 
consistent with cosmology out to values of the universal soft supersymmetry-breaking mass 
parameters mo,mi/2 that were larger than at tan/3 < 20. As a result, the discovery of 
sparticles at the LHC could not be 'guaranteed' in the CMSSM at large tan/3, unlike the 



case when tan/3 < 20 |2^. Since the recent BNL measurement favours qualitatively 
values of tan/3 that are not small, as does the LEP Higgs 'signal', and since the b sj 
constraint also begins to bite at large tan/3 even for > 0, it is important to understand 
the interplay of all these constraints. 

We find good compatibility between all these constraints for tan/3 > 10. Even if one 
generously allows a 2-cr downward fluctuation in the E821 discrepancy, one finds interesting 
upper bounds on mo and mi/2 that effectively extend the previous 'guarantee' of CMSSM 
discovery at the LHC to large values of tan /3. However, no such 'guarantee' can be offered 
to a linear e^e~ collider (LC) with centre-of-mass energy below 1.2 TeV. We discuss the 
uncertainties in our analysis associated with AQ,mi, and rrit [|. 

As already mentioned, our analysis is based on the one- loop calculations of ITS. In 



relating the masses of the sparticles appearing in the loops to the basic CMSSM soft 
supersymmetry-breaking parameters mo,mi/2, we incorporate the one-loop corrections for 
charginos and neutralinos. We also incorporate the leading two-loop electroweak correc- 
tion factor (1 — (4a;/7r)ln(m/m^)) [^, where m is a sparticle mass. We set the trilinear soft 
supersymmetry-breaking parameter = as a default, but we also discuss the consequences 
of varying it over the range — 2mi/2 < ^ 2mi/2- The b and t quark masses enter in our 
m/j, RGE and relic annihilation calculations. We use as defaults mb{mb)^M ~ ^-25 GeV and 
the pole mass = 175 GeV, commenting later on the changes as these mass parameters 



vary over our allowed ranges ±0.25 p8| , |33| , ±5 GeV. 

In our subsequent discussion, we consider the 2-cr range 75 x 10"^'' > Sa^ > 11 x 10~^° to 
be allowed by the E821 measurement |jl|, with the l-cr range 59 x 10"^° > (5a^ > 27 x 10"^° 
preferred. We interpret 75 x lO^^*' as a hard upper limit on (5a^, but models yielding 6a^ < 
11 X 10~^° should perhaps not be completely excluded yet. We note that a large amount of 



^Several other papers on the supersymmetric interpretation of the BNL measurement have aheady ap- 
peared during the last few days, and we comment on them at the end of this paper. 
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extra data have already been taken by E821, and that the present uncertainty will soon be 
reduced, which might have dramatic consequences. 

The LEP lower limit on the mass of the Higgs boson is rrih > 113.5 GeV, and the possible 
signal corresponds to nih = llStg:? GeV [||. This lower limit applies in the CMSSM, because 
the ZZh coupling is unsuppressed relative to the Standard Model ZZH coupling, unlike 
in general mixing scenarios possible in the MSSM. In the following, we display the range 
113 GeV < nih < 117 GeV. Given the uncertainties in the Higgs mass calculations 



choices of the MSSM parameters that yield slightly lower values of rrih might be acceptable, 
whereas values larger than 117 Gev are certainly allowed if one discards the LEP 'signal'. 

For b S7, we allow parameter choices that, after including the theoretical errors due to 
the scale and model dependences, may fall within the 95% confidence level range 2.33x10^^ < 
B{b — * s < 7) < 4.15 X 10^^. For the cosmological relic density, we allow the range 
0.1 < ^x^"^ — O-^- the lower bound is optional, as there may be other sources of dark 
matter, but the upper bound cannot be relaxed significantly. 

We display in Fig. the (mo,mi/2) planes for representative choices of tan/?, assuming 
yU > and Aq = 0. The regions allowed by the E821 measurement of at the 2-a level 
are (pink) shaded with solid black line boundaries. Also shown as black dashed lines are the 
regions favoured by at the l-cr level. We display the rrih = 113, 117 GeV mass contours as 
(red) dash-dotted lines, the dark (green) shaded regions are excluded by 6 — i> 57, the darker 
(red) shaded regions are excluded because the lightest supersymmetric particle is the lighter 
f, and the light (turquoise) regions are where 0.1 < n^h^ < 0.3. In panel (a) for tan/5 = 10, 
where there is no relevant constraint from 6 — * 57, we also show as a dashed line the lower 
limit m^± > 104 GeV. This excludes the tail of the cosmological region at large mo and 
small mi/2, where there is rapid annihilation through the h pole. For clarity, this chargino 
mass contour is not shown in the other panels, but its effect is similar |. 

We observe that there is remarkable consistency between the constraints from a^, rrih, 

b ^ s'-f and cosmology for tan (3 > 10, as also seen in panels (b, c, d) of Fig. |I[ Even given the 

uncertainties in the calculation of ruh, it is difficult to maintain consistency between and 

rrih for smaller values of tan f3. When tan (3 ~ 10, the other constraints are consistent with 

cosmology only for rriQ ~ 100 GeV, increasing gradually to rriQ ~ 170 GeV for tan/3 ~ 30, 

as seen in panel (b) of Fig. |l|. The favoured range of mo increases further as tan f3 increases, 

as a result of the cosmological constraint and the appearance, in particular, of the rapid 

XX ^ A,H annihilation process visible in panels (c, d) for tan (3 = 50,55, respectively, 

^We do not show the slepton mass constraint from LEP, which is weaker than the Higgs constraint in the 
CMSSM, and is weaker than the upper hmit 5an < 75 x 10^^°. 
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mi/2 (GeV) ^1/2 (GeV) 




Figure 1: The (mi/2,mo) planes for /i > and tan/3 = (a) 10, (h) 30, (c) 50 and (d) 55, 
found assuming Aq = 0,mt = 175 GeV and mb{mb)g'j^j = 4.25 GeV. The near-vertical (red) 
dot-dashed lines are the contours rrih = 113,117 GeV, and the near-vertical (black) dashed 
line in panel (a) is the contour m^± = 104 GeV. The medium (dark green) shaded regions 
are excluded by b s'-f. The light (turquoise) shaded areas are the cosmologically preferred 
regions with 0.1 < f^^^^ — -^^ dark (brick red) shaded regions, the LSP is the charged 
fi, so this region is excluded. The regions allowed by the E821 measurement of ttfj, at the 2-a 
level are shaded (pink) and bounded by solid black lines, with dashed lines indicating the 1-a 
ranges. 
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which allows mo < 800 GeV. The allowed ranges of mi/2 also increase as tan/3 increases to 
~ 50, where mi/2 ^ 900 GeV is allowed, falling slightly when tan/3 = 55 because of the 
rapid XX ^ A,H annihilation 

The allowed ranges of mi/2 are, however, much restricted if one uses the l-cr range for 
a^, with the maximum value being < 500 GeV. Indeed, combining all constraints and the 
l-cr range for a^, we find quite small allowed regions of the (mi/2,'n^o) plane centred on: 
~ (250,100) GeV for tan/? = 10, ~ (350,170) GeV for tan/? = 30, ~ (400,350) GeV for 
tan/? = 50, and ~ (400,500) GeV for tan/3 = 55 0. Typical sparticle masses corresponding 
to these choices are given in the Table. Comparing with the CMSSM physics reach for Run 
II of the Fermilab Tevatron collider, we see that the trilepton signature may be visible over 
some fraction of the allowed region of the (mi/2,'n^o) plane for tan/3 = 10, but not for the 
larger values of tan /3 studied here and in |3H . 
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m^± 






mil 


rrig 


m-g 


rrih 


rriA 
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250 
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99 
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145 


270 
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735 


790 


113 


475 
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350 
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315 


240 


385 
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875 


895 


114 


460 


40 


55 


400 


500 


170 


315 


315 


525 


665 


940 


895 


114 


450 


34 



Table 1: Typical sparticle masses for points in the (mi/2,'n^o) plane consistent with the 1-a 
range in and other phenomenological and cosmological constraints. For f, e, and t, the 
mass corresponds to the lightest eigenstate, which is mostly right-handed, and rriq corresponds 
to an average slight squark mass. 



As for the LHC, we have shown earlier that, in the absence of the LEP Higgs 
'signal' and the E821 value of a^, cosmology would allow mi/2 ^ 1400 GeV for tan/3 < 20 
in the coannihilation region. Studies by members of the CMS Collaboration have shown 
that at least some sparticles would be detectable at the LHC throughout this cosmological 
region More recently, however, it has been shown that the maximum value of mi/2 
increases to 1700 (2200) GeV for tan/? = 30(50), and that the 'funnel' of parameters allowed 
by rapid xX —^A,H annihilation also extends out to large mo and mi/2- These extensions of 
the CMSSM parameter space allowed by cosmology raised the spectre that the LHC might 
miss supersymmetry. 



•^Generically, we do not find consistent electroweak vacua for significantly larger choices of tan^, with 
our default values of mt , to;, and Aq . 

^Note that, for tan/3 = 10, we must relax the LEP Higgs constraint to remain compatible with the 1-a 
range for a^. 
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This is no longer a concern if the E821 lower limit 5a^ > 11 x 10^^° is confirmed. Fig. ^ 
shows the upper limits on mi/2 obtained as functions of tan (3 by combining cosmology with 
E821 or with the upper limit rrih < 117 GeV suggested by the possible LEP Higgs 'signal'. 
The constraint is somewhat stronger, but either would bring supersymmetry back within 
the range of the LHC. We also show in Fig. ^ the upper limits on itlq imposed by cosmology 
alone and in association with the mh or constraints. The rapid rise in the upper limit to 
mo from cosmology at large tan (3 is due to the appearance of the annihilation poles seen in 
panels (c) and (d) of Fig. |l|. Below tan/3 ~ 10, no independent limit on uiq is provided by 
the upper bound mh < llTGeV. For tan/? > 10, the limit on tuq is strengthened gradually 
as the rrih = 117 GeV contour slides down the coannihilation region, until this is offset by 
the shift in the cosmological region to higher mo as tan (3 is further increased. Eventually, for 
very large tan (3 > 45, the annihilation poles again allow very large values of mo. In contrast, 
the lower limit from always imposes a significant upper bound on mo. We conclude that 
the LHC will find supersymmetry, if the CMSSM is correct and the E821 lower limit holds 
up. 

It has been commented previously |37] that although discovery of the CMSSM could be 
'guaranteed' at the LHC if tan f3 < 20, there was no such 'guarantee' for a first-generation 
linear e'^e~ collider such as TESLA or the NLC with a centre-of-mass energy below 1.25 TeV, 
because of the extension of the cosmologically allowed region by coannihilation Q Since the 
E821 lower limit (5a^ > 11 x 10^^° excludes the 'tail' of the coannihilation region, the concern 
that such a first-generation linear e^e~ collider might miss supersymmetry is diminished. At 
the boundary of the region allowed by E821, the lightest detectable supersymmetric particle 
is the lighter stau fi. For our default choices mb{mb)s^M = 4.25, m^ = 175 GeV and Aq = as 
in Fig. I, and the test values tan /3 = 10, 30, 50, 55, we find that m^^ < 190, 320, 420, 580 GeV. 
We therefore conclude that a first-generation linear e~^e~ collider with centre-of-mass energy 
above 1.2 TeV would be 'guaranteed' to find supersymmetry within our CMSSM framework. 
A machine with centre-of-mass energy above 800 GeV would be similarly 'guaranteed' to 
find supersymmetry if tan (3 < 45. 

Fig. ^ for tan/? = 50,yU > and Aq = illustrates the effects of varying mb{mb)g^^^j 
between (a) 4.0 and (b) 4.5 GeV, and of varying rrit between (c) 170 and (d) 180 GeV. As 
one would expect, the changes in the Sa^ constraint are minor, being essentially associated 
with changes in the RGE and vacuum analysis. Also as expected, the m^ contours and the 
b —>■ S'-f constraint are rather similar in panels (a) and (b): the important change as one 
varies m^, is in the cosmological constraint. In particular, the rapid XX A,H annihilation 



'"'There was also no such 'guarantee' in the focus-point region p8|. 
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tan |3 

Figure 2: Upper limits on mi/2 and rriQ obtained as functions o/tan/3 for fi > 0, assuming 
mh{mf,)g^j^f = 4.25, rrit = 175 GeV and Aq = 0. We show the upper limits on mi/2 obtained 
by combining cosmology with the LEP Higgs 'signal' and the E821 lower limit on 5a ^, and 
the upper limits on rriQ imposed by cosmology alone and in association with either or the 
LEP Higgs 'signal'. 



'funnel' moves to lower mi/2 as m;, increases, reducing the combined upper limit on mi/2 

when mb{mb)s^M — 4-5 GeV, and increasing the upper limit on mo- However, our overall 

conclusions on the observability of the CMSSM at different colliders are unchanged. As 

seen in panels (c) and (d), the main effects of varying rrit are to move the m/j contours 

and the allowed cosmological region 0. As a result, the lower bound on tan/5 is relaxed for 

rrit = 180 GeV. However, the effects on the bounds on mi/2 and rriQ in Fig. ^ are again 

relatively minor. We do not display the effects of varying —2 x mi/2 ^ Aq < 2 x mi/2'- the 

main changes are in the allowed cosmological region, whose sensitivity to input assumptions 

were commented on previously ||28|, but the effects on the bounds on mi/2 and ttlq in Fig. ^ 

^Note also the black region in panel (c) of Fig. which is where we find no consistent electroweak vacuum. 
There are similar but smaller regions for larger m^, that are not shown. The size of this forbidden region is 
quite sensitive to the treatment of mj, a topic we leave for another occasion. 
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are again not very important, though the lower bound on tan/? may again be relaxed. 

In conclusion: we have combined the E821 constraint with other constraints on the 
CMSSM, including the LEP 'signal' for the Higgs boson, b ^ s'y and the favoured range of 
the cosmological rehc density fl^h"^. We find a high degree of consistency for tan/? > 10, 
and interesting upper bounds on mi/2 and rriQ. There is a corner of parameter space where 
the Fermilab Tevatron collider may find supersymmetry. On the other hand, discovery of 
supersymmetry is 'guaranteed' at the LHC, within our stated theoretical assumptions. The 
E821 6a^ constraint increases the chance that a first-generation e~^e~ linear collider will 
find supersymmetry, though there is no 'guarantee' unless its centre-of-mass energy exceeds 
1.2 GeV. 

The E821 measurement of 5a^ provides an important constraint on the CMSSM, and may 
already he the most promising positive evidence for it. With the prospects of a significant 
reduction in the E821 error bar in the near future, we may be living in exciting times for 
supersymmetry. 



Finally, for completeness, we comment on recent papers related to ours. 



Ref. 1 39 1 analyzes the BNL measurement but not the other constraints discussed here. It is 
suggested that sparticles may be produced at the Tevatron in Run II, but their observability 
is not discussed. In this regard, as discussed above, we are not very encouraged by previous 
studies P3[] . 

The most complete previous supersymmetric interpretation of the BNL measurement is 
given in |^ , which also includes some discussions of the LEP Higgs, b s'-f and cosmological 
constraints. The main difference between that work and ours is in the cosmology: the regions 
of CMSSM parameter space given in do not include all the coannihilation region, that 
extends in our calculations [H] up to mi/2 = 1400(2200) GeV for 0.1 < n^h"^ < 0.3 for 
tan/? = 10(50). Conversely, the 'focus point' region, that is now disallowed by the BNL 
and other constraints, is beyond the domains of the (wii/2,'"^o) plane that we plot. Our 
analysis for tan /5 = 50 also differs in the treatment of the direct-channel A, H poles, that 



are responsible for a large allowed region for ^^^/i in the analysis of |40|, but lead to narrow 



funnels in our analysis p8 



The prospects for dark matter detection in the light of the BNL measurement are dis- 
cussed in [^, where LEP, 6 — > 57 and fi^^/i^ are also taken into account. The current LEP 
constraint > 113.5 GeV was not used. This applies in the generic MSSM for tan/? < 10 
and in the CMSSM at essentially all tan/5, excludes a large range of mi/2 and requires 

We are not in a position to compare treatments of 6 ^ sj. The 



> 150 GeV 128 
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mi/2 (GeV) nii/2 (GeV) 

Figure 3: Comparison between the (mi/2,'n^o) planes for tan/? = 50, /i > and Aq = 0, 
with different values of other input parameters. Panels (a) and (b) are for rrit = 175 GeV, 
'fnbijnbj^M = 4.0 andi.d GeV, respectively. Panels (a) and (d) are for mb{mb)'g^^ = 4.25 GeV 
and nit = 170 and 180 GeV, respectively. 
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analysis of Q^h"^ in PT| is based on DarkSUSY which does not include all the effects at 



large tan (3 that we discussed in |Q and here. 

We use the same formulae as to implement the BNL constraints on mo and mi/2, 
with which we agree quite closely. Ref. also discusses the LEP Higgs 'signal' and makes 
qualitative comments on supersymmetric dark matter, but does not discuss b sj. 



Ref. discusses implications for the unconstrained MSSM and for gauge-mediated 
models of supersymmetry breaking. 



Ref. 1^ discusses the supersymmetric interpretation of in relation to models of 
neutrino masses and the observability of /i ^ 67 decay. 

We 



We note that anomaly- mediated supersymmetry breaking is considered in [40, 43 



also note that non-supersymmetric interpretations of the E821 result are discussed in |46 
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